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1. INTRODUCTION 

Real-life applications investigate the meaning of fuzzy as a generalization of crisp common sense 
because it is a solid instrument for modeling the vagueness, and in specific, to treat uncertainty with a 
mathematical model [1]. In many real well-determined dynamic issues, a system of ordinary or partial 
differential equations may represent the mathematical model. On the contrary, fuzzy differential equations 
(FDEs) are a valuable tool to model a dynamic system that is ambiguous in its existence and comportments. 
Since the FDEs have been used widely to model scientific and engineering problems, they have become a 
popular topic among researchers [2]. There are many practical problems with the solution of FDEs that 
satisfy initial [3] or boundary [4] values conditions. The main reason why finding the approximate solutions 
to the fuzzy problems becomes necessary is that most of the problems are too complicated to be solved 
exactly, or there are no analytical solutions at all [4]. Hence, FDEs will be suitable mathematical models for 
dynamic systems where complexity and ambiguity occur. For this reason, we may find FDEs exist in several 
fields of mathematics and science, including population models [5]-[7] and mathematical biology and 
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physics [8]-[10]. As an alternative to the analytical solutions to such problems, approximate solutions such 
as homotopy perturbation method (HPM) and variational iteration method (VIM), and Adomian 
decomposition method (ADM) are listed as some of the approximate analytical methods [11]. 

In the last decade, some researchers have started to explore the numerical solutions for two-point 
fuzzy boundary value problem (TPFBVP) [12], [13]. Semi-analytical approaches have been used over recent 
years to overcome linear TPFBVP by various methods [14]—[16]. He implemented the HPM in 1999, and the 
method was applied to a wide range of mathematical and physical problems [17], [18]. This method provides 
the solution into components of a short convergence series, which are elegantly determined. HPM is now 
known as a standard tool for overcoming all kinds of linear and nonlinear equations, such as differential or 
integral equations. Another significant advantage is that the measurement size can be reduced while 
increasing the exactness of approximate solutions, so it is regarded as a robust process. Along with the HPM, 
VIM was introduced and also proposed by He [19]. This approach differs from specific classical techniques 
utilizing which nonlinear equations are quickly and accurately resolved. VIM has been used in many physics 
and engineering sectors recently [20], [21]. This approach is helpful for directly solving linear and nonlinear 
problems with n“ order boundary value problems (BVPs) without reducing them to a BVP system. It has 
been reported by many authors, such as [22], that VIM is more robust than other analytical approaches, like 
ADM and HPM. Compared to HPM and ADM, where computer algorithms are commonly used for nonlinear 
terms, VIM is used explicitly without any nonlinear terms requirement or restrictive assumptions [23]. 
Without restrictive assumptions, the VIM solves differential equations that can change the structure of 
solutions. In VIM, the calculation is simple and straightforward [24]. The VIM overcomes the difficulty of 
measuring Adomian polynomials [23], a significant advantage over ADM. 

According to [24], one significant disadvantage of VIM is that the terms obtained are longer than 
those obtained by decomposition and perturbation methods. For this reason, we are seeking and investigating 
the proposed approximate analytical solutions for nonlinear FDEs by HPM and VIM for comparison 
purposes. This study will develop an innovative approach to modifying the nonlinear TPFBVP based on the 
framework of fuzzy problems. This modification is tested on two existing FDEs and compared with the exact 
solution and the numerical solution. A comparative study will be given to show the capabilities of the 
proposed methods. According to the results, the modified schemes were found to be feasible and more 
accurate. 


2. PRELIMINARIES 

In this section, we provide some fundamental concepts and definitions that are necessary for this 
work. This includes some propositions, properties, and explanations of fuzzy sets and numbers and FDEs that 
will be used later in this work. 
Definition 2.1 [25]: the relation: 


0, ifx<a 
x-a ae ee 
Baa’ ifa<x<B 
ux; a, B,y) = y-x fg 
y iff <x< 
y-eP 7 
0, ifx >y 


is a form of the membership function for a trapezoidal fuzzy number u = (x; æ, ß,y), which is presented as 
shown in Figure 1. 


> 
X 


Figure 1. Triangular fuzzy number 
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And its r-level is: [u], = [a+ r(B-—a),y — r(y-ßp)ļ} re [0,1]. 
Definition 2.2 [25]: The relation 


0, ifx<a 
x-a . 
Boe" ifa<x<f 
u(x; a, B,y, 5) = 41, if B<x<y 
6-x Fi 
zy’ ify <x <ô 
0, ifx >ô 


Is a form of the membership function for a trapezoidal fuzzy number u = (x; a, ß,y, ô), which is presented 
as in Figure 2. 


Figure 2. Trapezoidal fuzzy number 


This can be used to describe an r-level set of the trapezoidal fugitive number as; 
lAl = [$ -a)r + æ, — (8 — y)r] 


This paper describes the class of all fuzzy subsets of R is being marked by Ë that satisfies the characteristics 
fuzzy number properties [26]. 

Definition 2.3 [27]: Let f: R > E, f(x) is called fuzzy function if E is a set of fuzzy numbers. 

Definition 2.4 [28]: The r-level set defined as [f(x)], = [Fr Fs r)| x € K,r € [0,1] for a fuzzy 
function f:T > E where T c Ë. A fuzzy number is more effective than the r-level sets as representational 
types of fuzzy sets. Fuzzy sets can be described based on the resolution identity theorem by the families in 
their r-level sets. 

Definition 2.5 [29]: If f: X > Y is function induces another function f: F(X) > F(Y) For each interval, U in 
X is defined by: 


Supyer-1(4) U(x), if y € range (f) 


rwo- „if y ¢ range (f) 


This is recognized as the theory of Zadeh extension principle. 

Definition 2.6 [30]: Let Z = X © ¥ be the H-difference of the fuzzy numbers % and ¥ if the fuzzy number Ž 
exist with the property% = y + Z. 

Definition 2.7 [13]: If f: I > E and ty € I , where I € [a,b]. f’ is said to be Hukuhara differentiable at to, if 
there exists an element [F], € E to be small enough for all h > 0 (near to 0), exists f(tp + h;r) © 


f (tosr), FC(to;r) © F (to — h; r) and limits in metric (Ë , D) are taken and exist in such a way that 


fatin OH tot) = Jim Fton OR torte) 


f (to) = lim 
f ( 0) h>0+ h>0+ 


For more details, see [30]. 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1916-1928 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1919 


Definition 2.8 [31]: Settle f : I > Ë and ty E I , for € [a,b]. f™ is said to be Hukuhara differentiable 
x € Ë, if there exists an element [f w € Ë to be small enough for all h > O(near to 0), exists 
FTVM(ty + hir) O FOD (ta r), FOD (toir) © fY (to — h;r) and limits in metric (Ë, D) are taken 
and exist in such a way that 


FO) = lim FEV (tothirjOFfeM(toir) _ lim FEY (torIOFOYV(to-hir) 
p h>0+ h h>0+ h 


There is a second order of the derivatives of Hukuhara for n = 2 and equivalent to f™. 
Theorem 2.1 [31]: Let f: [t + a,b] > E be Hukuhara differentiable and denote 


Fol, Foro] = [enren] 


Then we can define the differentiable boundary functions f'") and f'(t; r) can be written in the n™ order of 


n" fuzzy derivatives. 


[Feo], = G (t; a )h vr € [0,1] 


3. DESCRIPTION OF THE FUZZY HPM 
The overall HPM structure for solving crisp nonlinear TFBVP is mentioned in [17], [18]. Consider 
the defuzzification of the following general n* order TPFBVP [15]. 


IMO) = f (650.7 OHO, FOO) + GO, tE [toT] (1) 


te = ZO, F (to) = ay, oo I(t) = a, (2) 


HT) = FO, FT) = BO, ...., GODT) = POD, 
To solve (1) by using HPM, we need to fuzzify HPM and then defuzzify it back to (1) as in [15]. According 
to [32], the HPM and for all r € [0,1]. HPM must ensure the convergence of the HPM solution series 


function through the correct choice of initial guess and the auxiliary linear operator [33]. From [15], HPM 
form for solving (1) is given by solving the lower bound as follows: 


Ly [yo Ct; r)— yo(t; r; weiG@)| = 0, 
p°:4 ltor) = [æo] y (torr) = [aa]. (torr) = [ana] 
y(Tir) = [Bo] oy’ ir) = [Bi] sy? Tir) = [Brea]. 


Ly [ya (ts) + volt: r; Dt cs(r))| — F (t, (t; r, Zhi csCr))) 
—g(tir) = 0, 
Yı (tos) = 0, y1 (tor) =0, es P Cto; r) =0, 
y (Tir) = 0, y (T; r) = 0, .., y1 (T; r) = 0. 


Lny2(t;r) -F (t, 0 (t:r, £ c())) = 0, 
p”: Y2(to;r) = 0,,¥2'(toir) = 0, ma eget) =0, 
y2(T;r) = 0592'S Oe yo’ OE) = 0. 


LyYnes(tit) — F (t,0(ti7, Dh es(7))) = 0, 
pers Yn+i(tos T) ot 07, Yn+1 (tor) =0, no Yna P (Cto; r) =0, 
Yne1(T37) = Or Yni (Tir) = 0, aes Yn (Tir) = 0. 
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Similarly for the upper bound 


Laly r) -Yol r; D1 e5(r))] = 0, 


0:2 H(tosr) = [Bol F (torr) = [Eales We? (tor) = [Rail 


p : 


ensko Cosas Oy) = Baes 


La [yy Er) + 7t r; ELE] - G (t, olt; r, Eia E), 
1. —g(tir) = 0, ss 
y, (oir) = 0,y¥, (to; r) = 0 , a Vy i (to; r) = 0, 
7, (Tr) = 0,7, (Tr) = 0,9, * OC) = 0. 


LREN — G (ihr, Dis -(r))) = 0, 


p : Y, (to; r) =, 0,¥, (toir) = 0 , gn tae) = 0, 
F(T; r) = 0,7, (T;r) = 0,5, P (T; r) = 0. 


A T A O 
Faaa ltor) = 0 Fpp (toir) = 05 Fpa oir) = 0, 
Faaa Tr) = O Fpa TIT) = 0, a Pp PT; r) = 0, 


n+1, 
p : 


where V(t) = H(t), 9'(t), Y” (©, ... 9" P(t) and the initials guessing F(t; r) is given in [16] for all r € 


[0,1] and then the approximate solution is given by setting p = 1 as in (3): 
ICE r; Loar Cs(7)) = Sin (ti r; Dr ĉs(r)) = LiPo Vile r; Leas &Cr)) 
Therefore, the exact solution of (1) can now be obtained by setting p=1: 


Y(t;r) = lim F(t r; Yea ĉr)) = lim{Qi=o pVAGT; Xs-1 č r))} 
= Vico Vi(ts r; Deas č&(r)) 


4. DESCRIPTION OF THE FUZZY VIM 


(3) 


(4) 


VIM general structure for resolving problems TFPBVP is stated in [19]. In order for the (1) to be 
solved using VIM, we must fuzzify and defuzzify the VIM as defined in (2) [34]. According to VIM in [22] 


and for all r € [0,1] we rewrite (1) in the following correction functional forms: 
Wea Glee) Hy (61 Tere) + 


KAn) feon + F (7 a) + gain) dn, 


Yip; r; Ès=1 Cs (r)) 


=F r ELT) + f, ACs) ip (nr) + G (n? (n; n) +g(n; n| dn, 


(6) 


where i= 1,2,..., r € [0,1]. The Lagrange multiplier is A(t;7) which can be optimally defined via 


variational theory [33]. Now we let 


R= P(n? onr)) + g(n;r), 


GH, =G (n? onr)) +g(n;r), 
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where F, and G, are nonlinear operators including the nonlinear terms F and G and the inhomogeneous term 
g(;7r). In the following forms we will rewrite (5) and (6) as: 


Viva (t; r; Elaa c (r)) = yu ters Da os (7)) + Sp Actin) {i Cir) + RH} an, 
Vi (tir Doha G(r) =H, (ti r; DE Ger) + So Acti nd FO Gir) + GH} an, 


where restricted variation is F, ie ô F = 0 [21]. The general multiplier A(t; n) applied to (1) according to 
[23] can be described in the following: 


Syin EEG) = Sylt r Ei e(r)) + 5 Alt m {y r) + E F} an 
= y (Grist e(r)) = öyltr DIP cs) +6 f Altin) {yi 0; r)} an. 


According to [23], we obtain the followings by integrating by part: 


Yili r Esr Ge) = [1 -40 ]Sy (tir Tee) + Ee? Syt; r) + 
fy Alt; MP dy; (t; r; Z2 cs (r)) dn. 


The following stationary conditions can therefore be reached [35]. 


am™ = 0, 
1-400 =0, 
A(t) = 0,k = 0,1,2, ..,n— 2. 


Under these terms and the order of the (1) the general Lagrangian multiplier can be calculated as [35]: 


—1\nfín-pn-1 
atm = Pe) 


(n-1)! 


Therefore, all the above parameters A(t; n) and Folt; TEs (r)) will be easily obtained in the series of 
approximations of VIM. The exact solution can therefore be obtained: 


F(t r) = lim y(t; r; Leas C5 (7). 


To determine the č (r), we use the same HPM technique as in [15] by substituting these constants with the 
initial estimates in series solution function and then using the boundary conditions of (1) to determine the 
values of these constants for each fuzzy r-level set. 


5. NUMERICAL EXAMPLES 

The approximate solution by HPM and VIM in the next examples are obtained by formulating the 
given equations as presented in sections 3 and 4. The formula is then solved and analyzed by using 
Mathematica 11: 
Example 5.1: Let us consider the following nonlinear second order TPFB VP: 


F'OPO =t*+2, 90) =4 590) =b (7) 


where @ and b are triangular fuzzy numbers having r-level sets [0.17 —0.1,0.1—0.1r] and 
x 2) 
[0.9 + 0.1r,1.1 — 0.1r] for allr € [0,1]. The linear operator in compliance with section 2 is L} = i with 


the inverse operator ee and the initial guesses for all r € [0,1] are given by: 


yoltir) =alr)+ o (r)t, 
a (tr) =¢C,(r) + e,(r)t. 
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5.1. HPM formulation 
According to section 3, the values of & (r) = [0.17 — 0.1,0.1 — 0.17] and the homotopy functions 
of (7) are 


H(t,pir) = (1 -PL [yn — yo (tirica@))] +p ave n+ (yen) -e+ 2| = 0, 
H(t,p;r) = (1 -p)L |y; r) — yo (; r; (r) )| +p vee r)+ (vee r) — (tê + 2| = 0, 
H(t,pir) = (1- p)L[ytir) -¥,(tsr;&))] + p[ Layer) + (Wer) -e +2)] =0. 


The n components are specified as the HPM in section 3 of ¥,(t;r) for k =1,2,..,n and r € [0,1] 
calculated by evaluating the lower limit as follows: 


p’: {yo (ts r) = 0.1r — 0.1 + cœ (r)t, 

pis y(t;r) = L, (-» (t; rie) +t + 2) 
y0; r) = 0,y, (ir) = 0. 

pe a r) =—-L,* Ba (t; rio (r)) yo te riel) )| , 


y2(0;r) = 0, y2(1;r) = 0. 


Ja pane r) = -L Eii ye (Eiri (7) Yni (HT EOC), 


(8) 
Yr(0; r) = 0, Yk (1; r) = 0. 
For the upper limit, the values are obtained in the same way as follows: 
p°: Ct; r)=0.1-01r+¢,(r)t, 
ane ar S 2 
ET yr) =L i |- a(t; T; C2 (r))) +t*+ 2], 
y,(0;r) = 0,y, (sr) = 0. 
ee P 
y,(0;r) = 0,y,(1;r) = 0. 
= > -1 aq — = p _ 
waa [Yr Er) = -L2 tec Y(t r; C2 (r)) Yn-1-k (t; r; C2 (r))], 
po 2 = (9) 
Yr) = 0, y, (1; r) = 0. 
Evaluating (8) and (9) to obtain fifth order HPM series solution in the following form such that 
Š: (t; r; G(r) = Polti r, G(r) + Xi- Hi(67 &lr)) = Hore) (10) 


Now to obtain the values of ĉ (r) for allr € [0,1] , we solve the nonlinear series solution of (7) from the 
boundary condition [0.9 + 0.1r,1.1 — 0.1r] then we substitute the values of ĉ (r) again in (10) to obtain 
fifth order HPM series solution. Since (7) is considered without exact analytical solution, to show the 
accuracy of fifth order HPM approximate series solution §,(t; r) for all r € [0,1], the residual error must be 
specified: 


[E(t)], = Sst r O + S(t; r) E) -tt + 2, 
Then the fifth order HPM series solution is presented in the Tables 1-2 and Figure 3. According to Tables 1 
and 2 and Figure 3, we concluded that the fifth order HPM approximate solutions of (7) for all t € [0,1] 


and r € [0,1] fulfill the patterns of fuzzy numbers in the form of a triangular fuzzy number. 
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Table 1. Approximate solution y(t; r) of fifth order HPM at t = 0.5 andr € [0,1] 


r &(r) y(0.5;r) [E (0.5)], 

0 —0.012233517979781578 0.14366068333554236 4.183283991154862 x 1077 
0.25 —0.010209834659421881 0.16997259607001605 2.946362775008149 x 1077 
0.5 —0.007505097713710612 0.19646322624295728 2.485028610110795 x 1077 
0.75 —0.004110208069923338 0.22313542853095245 2.726826934945636 x 1077 

1 —0.000016001865820840 0.24999204985486695 6.342108563664461 x 1077 


Table 2. Approximate solution y(t; r) of fifth order HPM at t = 0.5 andr € [0,1] 


r &(r) y(0.5;r) [E(0.5)], 

0 0.023538969063634813 0.35931908435924000 0.0000211620477442997 
0.25 0.016554919244886100 0.33169666727607905 0.0000108129304682291 
0.5 0.010307281584441760 0.30426986657956345 0.0000048897414700183 
0.75 0.004786740073139446 0.27703592414404127 0.0000018877716657511 

1 —0.000016001865820840 0.24999204985486695 6.342108563664461 x 107” 

, __5-Order HPM Solution 
1.0 
0.8 
2 
B 06 
T 
> 
® 
~ 0.4 
L 
0.2 
0.0 4 


0.25 0.30 


y(0.5;r) 


“0.35 


Figure 3. HPM approximate solution of (12) at t = 0.5 and r € [0,1] 


5.2. VIM formulation 
The variational formula of this problem is given in accordance with section 4 as (11): 


+ f Aen) {yi (n; r;c(r)) + Yi (n; no) == 2} dn, 


+ falt; n) D'r; T2 (r)) + y;(n; r; 2 m) -—n*- 2} dn. 


Viet G r; c(r)) =y; (t; r; c (r)) 


Var (tir; c(r)) = y;(t; r; €2(r)) 


The fourth-order VIM series solution is obtained in the form 12: 


Siete) = Eto Filtr &r)) =J; r) 


(1) 


(12) 


The Lagrangian multiplier of (11) is described in section 4 such that A(t; n) = n — t. Now in order to obtain 
the values of €,(r) for allr € [0,1] , we solve the nonlinear series solution of (7) from the boundary 
condition [0.9 + 0.1r, 1.1 — 0.17] then we substitute the values of ĉ (r) again in (12) to obtain fourth-order 
VIM series solution. For (1), the following residual error is described in order to demonstrate VIM accuracy 


in approximate fourth-order solution without an exact analytical solution: 


EO = |S." (ter) + ent - 2. | 


In Table 3 and 4 and Figure 4, the fourth order VIM series solution is presented: 
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Table 3. Approximate solution y(t; r) of fourth order VIM at t = 0.5 and r € [0,1] 


r &(r) y(0.5;r) [£(0.5)], 

0 —0.012227411850487582 0.14366373723600040 1.268180005453700 x 10717 
0.25 —0.01020310975429537 0.16997595759445022 3.399442879969960 x 10711 
0.5 —0.007496876529468979 0.19646732842074352 6.365820370834996 x 10711 
0.75 —0.004099110842952412 0.22314095416570680 8.448242105885129 x 10-1 

1 5.9121871734293 x 1078 0.25000002937629756 2.485910088889653 x 1071? 


Table 4. Approximate solution y(t; r) of fourth order VIM at t = 0.5 and r € [0,1] 


r c(r) y(0.5;r) [E(0.5)] 

0 0.023538969063634813 0.35935808190205626 7.0343751656931630 x 107? 
0.25 0.016554919244886100 0.33172332776916430 3.9969035997566850 x 10°? 
0.5 0.010307281584441760 0.30428777924960454 2.0336278494514910 x 107? 
0.75 0.004810786651353893 0.27704784239391295 8.636787535154511 x 1071? 

1 5.9121871734293 x 1078 0.25000002937629756 2.485910088889653 x 1071? 

, __4-Order VIM Solution 
1.0 
0.8 
£ 
® 
N 0.6 
T 
> 
© 
— 0.4 
A 
0.2} 
0 of 1 1 à š J 
0.15 0.20 0.25 0.30 0.35 
y(0.5;r) 


Figure 4. VIM approximate solution of (13) at t = 0.5 and r € [0,1] 


According to Tables 3-4 and Figure 4, we concluded that the fourth-order VIM approximate 
solutions of (7) for all t € [0,1] andr € [0,1] comply with the fuzzy numbers of properties as triangular 


fuzzy number. 


Example 5.2 [35]: Consider this non-homogenous second-order non-linear TPFB VP: 


P'O ++I = VO) + FO, € [0,1] 


F0; r) = ã& 90) =b 


(13) 


where @ = l =) + iE) + ] and b = E =) + = =(=*) + = for allr € [0,1]. According to 


’4\ 10 


[30] the fuzzy function f(t) have the following defuzzification: 


aa aa a GEIE 


Then the corresponding analytic solution of (13) is given by 


eagar a 


by following the initial guesses from Example 5.1, we have: 


r-2 1 


yoltir) =i) ttet 


7 tr) = TE) +2 40,(r)t. 
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5.3. HPM formulation 
From Section 4 the approximate solution of (13) can be determine by HPM from the followings: 


1 (r-2 


p°:{yo(t:7) =2(—) +++ a(t, 


4 


pi: yı(t;r) = Loe? pot r)i+ (v (t; r;o(r))) + f(t; |, 
yı(0;r) = 0,y,(1;r) = 0. 


(16) 
k-1 i 
pktt; Yr+1(t; r)= Lat yxlt; r) + > > Yk-1-i (t; r; c(r)) Yj (t; r; c (r)) Yi-j (t; r; c(r)) : 
i i=0 j=0 
yx (0;r) = 0,y,(1; 7) = 0. 
Similarly for the upper bound 
0. fF, (6:7) = (aa) +g? (r)t 
p a Yo ’;r = 4 10 4 C2 F , 
= S = = 3 = 
pt: y, (tir) = L [z,c r) + (a r; ca(r))) + f(t; r)| f T 


y,(0;r) = 0,7, (1;r) = 0. 


peri, Pen ED = LaF, (G7) + BIS jeo Peni (sroM)F (iren) 7 (trat) 
Y0; r) = 0,7, (131) = 0. 


Next, the determination of the values of ĉ (r) for allr € [0,1] can be done by solving the nonlinear seventh 


= -1 
order HPM series solution of (17) subject to the fuzzy boundary condition y(1;r) = (=) + = and 


-1 
y(1;r) = TE) + = of (13). Then, we substitute the values of ĉ (r) again in the seventh order series 


solution of (13) to obtain seventh order HPM series solution in Table 5. 


Table 5. Approximate solution of seventh order HPM of (13) att = 0.5 Vr € [0,1] 


r olr) T(r) y(0.5:r) y(0.5;r) 

0 —0.44985575506896760 —0.04891373031802654 0.0016323384080985234 0.30163258314132690 
0.25 —0.42489220652617230 —0.07406371025083800 0.0203823093927938160 0.28288260102272034 
0.5 —0.39991573435157324 —0.09920243697696327 0.0391322807152672800 0.26413260400295260 
0.75 —0.37492630650497016 —0.12432965360921329 0.0578822454883946360 0.24538262188434600 

1 —0.34992391523689814 —0.14944512128297516 0.0766322059390083700 0.22663263231515884 


5.4. VIM formulation 
The VIM to solve (13) is in compliance with section 4. The formulation shall be given as in (18): 


Vier (¢ r; olr)) =y; (& r; c (r)) 
y" mron) -y (tirer) 
-y (mro) - fr) 

Viale T; Ca (r)) = SAGs r; C2 (r)) 
„ina 0720 erao) 
Janam) -ftr) 


+ fin- t] dn, 


(18) 


By following VIM formulation and analysis in Example 5.2, third-order VIM series solution is given in the 
following Table 6: 
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Table 6. Approximate solution of third order VIM of (13) att = 0.5 Vr € [0,1] 


r olr) c (r) y(0.5;r) y(0.5;r) 

0 —0.4500000805284215 —0.05000128774459855  0.076632688439284070  0.30163199259313456 
0.25 —0.4250000216571309 —0.07500108858114996 0.057882686548291370  0.28288209737819590 
0.5 —0.3999999813340253 —0.10000090368612918  0.039132675010097500 = 0.26413219449035540 
0.75 —0.3749999595833066 —0.12500073372107165  0.020382653821232450 0.24538228361308428 

1 —0.3499999563806096 —0.15000057928658794 _0.001632623001234793 0.22663236446104182 


In order to show efficiency of VIM and HPM methods in solving (13), numerical comparisons of the 
accuracy generated by RKHS in [35] are presented in Table 7. These comparisons are conducted at t = 0.5 


and various r fuzzy level that belongs to [0,1] by computing the mean of the average error. Here, E(t;1) is 
the mean of the average error between E(t; r) and E(t; r) such that 


E(tr) = ær) -yer |, 


Elt; r) = |Y r) -7t r). (19) 


Table 7. Numerical comparison of mean errors E(t;r) att = 0.5 
T Third-order VIM Seventh order HPM RKHS [33] 
0 1.267224325474103 x 1077 2.397203478878614 x 1077 3.87851 x 1076 
0.25 1.690374040186348 x 1077 2.369393627074407 x 1077 3.31263 x 1077 


0.5 2.251369623368360 x 1077 
2.948819071355024 x 1077 
1 3.780944354508757 x 10-8 


2.521985155955197 x 1077 
2.508378323807635 x 1078 
2.630438642283513 x 10-8 


1.84790 x 1077 
1.34477 x 1077 
1.09977 x 1077 


According to Tables 6 and 7 and Figure 5, we concluded that the third order VIM and seventh order 
HPM will successfully provide the approximate solutions to (7) for all t € [0,1] and r € [0,1]. This is shown 
by the results that comply with the fuzzy numbers of properties as triangular fuzzy number. 


et rr T 


tof * + 


* VIM 
* HPM | 
* Exact J 


t-level sets 
w 
X- 


Figure 5. Exact solution of (13) compared with third order VIM and seventh order HPM 


approximate solution of (13) at t = 0.5 andr € [0,1] 


6. RESULTS COMPARISON 


In this section, we present a comparative explanation between the solutions of nonlinear TPFB VPs 
obtained by HPM and VIM as we illustrate from Examples 5.1 and 5.2: 


The initial approximation guesses in HPM, and VIM are obtained in the same way. 
— The construction of VIM formula to solve nonlinear TPFBVP is faster and easier than HPM because 
HPM takes the advantage of the small parameter p € [0,1] that makes HPM suffers from the cumbersome 


work needed for the derivation of for nonlinear terms. This will increase the computational work 
especially when the degree of nonlinearity increases. 
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— As mentioned in section 4, VIM is used directly without any requirement or restrictive assumptions that 
the nonlinear terms make the series solution longer and time consuming in CPU. The third order VIM 
solved (7) within 31.2264607 seconds for reach r-level values. In HPM the use of embedding parameter 
p is decomposed for the nonlinear terms making the series solution shorter than the solution of VIM with 
less time consuming in CPU. The seventh order HPM solved (7) within 1.4882989 seconds for all 
r-level values. For illustration, the nonlinear term y? in (7) which has been decomposed to 

n-i Ve Yn-1-k in HPM formulation but in VIM formula we substitute y. 

— From the results obtained by HPM in Tables 1-2 and VIM in Tables 3-4, we conclude that VIM provides 
a better and more accurate solution than HPM, with less order of series solution for both t € [0,1] and 
r E [0,1]. 

— Finally, both Figure 3 and Figure 4 show the solution of (7) by using HPM and VIM respectively for all 
t € [0,1] and r € [0,1] satisfy the fuzzy numbers properties in the form of triangular fuzzy number. 


7. CONCLUSION 

In this work, approximate analytical methods have been used for nonlinear TPFBVPs to achieve an 
approximate solution. Two schemes, HPM and VIM, were developed and reformulated to approximate the 
nonlinear TPFBVP solution. Numerical examples, including nonlinear TPFB VPs, demonstrate the efficacy of 
these approaches. Forgiven nonlinear TPFBVPs, a comparison of HPM and VIM results was presented. The 
comparison shows that VIM convergence is faster and provides an improved solution, particularly for the less 
approximate terms, nonlinear TPFBVPs over HPM. Even though these equations are without exact analytical 
solutions, the exactness of both HPM and VIM can be calculated from nonlinear TPFBVPs. The VIM 
Lagrangian multiplier for the nonlinear TPFBVPs is equivalent to the value for all the r-level sets. All 
outcomes of the experiments with HPM and VIM are achieved using a triangular shape to acquire the 
properties of the fuzzy numbers. 
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